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We formulate a freeze-out criterion for ultra-relativistic heavy ion collisions in terms of 
the pion escape probability from the collision region. We find that the increase in pion 
phase-space density from SPS to RHIC reported at this conference has a small influence 
t-h ■ on particle freeze-out because of the small tttt cross-section. Our treatment takes into 
account dynamical expansion, chemical composition and momentum of the escaping par- 
ticle for particle freeze-out. It supports a freeze-out at rather low temperature — below 
100 MeV — and earlier decoupling of high p± particles. 

o : 

A precise formulation of the condition under which particles decouple ("freeze-out") 
from a heavy ion collision is important for understanding the dynamical origin of the final 
state. The observation of STAR [jl] that the pion phase-space density increases at RHIC 
i— h \ indicates that freeze-out is not characterised solely by a universal value of pion phase-space 
density, in contrast to a recent suggestion [[|. Also, freeze-out is not characterised solely 
by a universal value of the total particle density since the volume from which particles 
> . decouple shows a non- monotonous dependence on collision energy ||. 

It has been argued 0] and demonstrated || that the chemical composition of the system 
must be taken into account. For pion freeze-out, the nucleon and anti-nucleon density is 
more important than that of pions since the pion-nucleon cross-section is larger than the 
pion-pion one. With increasing collision energy the pion abundance and their contribution 
to scattering grows, but it is argued that the momentum- averaged mean free path of a 
pion at freeze-out stays constant ||. 

This argumentation still ignores the influence of global dynamics on the freeze-out crite- 
rion. A strong expansion makes the particle density decrease fast and if the characteristic 
time scale for density decrease r exp = (d fJ- u fl )~ l is smaller than the mean free time between 
collisions T scatt a particle is likely to seize interacting, i.e., to freeze-out 0. 

In our work we study — in addition to these two aspects — how the probability for a 
particle to decouple depends on its momentum 0. We formulate freeze-out as condition 
for an individual particle: a particle freezes-out at the time of its last scattering. This 
is in contrast to the "fireball freeze-out" which is understood as a boundary in space- 
time between interacting matter and free-streaming particles. In the extreme case of 
a hydrodynamic model the "fireball freeze-out" is usually treated by the Cooper-Frye 
prescription []7| in which all particles decouple along the same three-dimensional hyper- 
surface when the medium fulfils certain condition. Our treatment goes beyond Cooper- 
Frye: particle freeze-out depends on its type, position, and momentum in addition to the 
characteristics of the medium. This corresponds to freeze-out as described by a cascade 
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event generator. 

The freeze-out criterion which we explore is based on the escape probability || 

V(x,p, t) = exp ^— J df 1Z(x + vf , p) \ , (1) 

which determines the probability that a pion emitted with momentum p from the position 
(x, t) escapes from the medium without further interaction. Here, the scattering rate 
lZ(x,p) is integrated along the trajectory of the pion. Freeze-out is assumed to occur 
when V(x,p, r) reaches a universal value of order one. 

The role of collective expansion: To illustrate how collective expansion affects particle 
freeze-out, we consider the case of particles with vanishing momentum in the centre of 
a fireball with longitudinally boost-invariant and transversely linear expansion velocity 
profile of the form 

= (cosh 77 cosh 774, cos0 sinh?7£, sin0 sinhr/t, sinh?7 cosh 774) . (2) 

Here, 77 = Arctanh(z/t) and r] t = £r where £ = 0.08 fm _1 is a typical value for the 
transverse expansion gradient |§. At the time r cm at which the pion is emitted, the 
density of scattering centres decreases by the rate 

idp 

p dr 

A power-law decrease of density with time is consistent with (£3) if 



dX = — + 2£- (3) 



T 



p{r) = Pem ( — ) , a = 1 + 2£r cm . (4) 



Under the assumption 1Z oc p we obtain for the opacity integral of a particle of zero 
momentum 

d T TZ(T) = —^-T em = , (5) 
a - 1 2£ 

where 7£ em is the scattering rate at r cm . Technically, the calculation of this opacity 
integral for particles of non-zero momentum p is more complicated, since one has to follow 
the propagation of the particle through layers of different density. For a quantitative 
statement, this requires a model of the space-time evolution of the fireball. Qualitatively, 
however, the general result will be consistent with the main feature of the case p = 
considered here: as seen from (|5|), freeze-out from a denser fireball (larger lZ em ) is possible 
if it is compensated by a stronger expansion gradient such that the value of opacity integral 
stays constant. 

The role of chemical composition: To determine the value of 7£ em for collisions at \fs = 
17AGeV (SPS) and y/s = 130AGeV (RHIC), we calculate the scattering rate 

^(p) = E /d 3 A^(fcK(s)k-^l, (6) 

i 

where a[ is the cross-section for pion scattering on species i. For the collinear cross-section 
we use the parametrisation 

«M = Ett. <*. ".I I*. (2Si I + l) ±; s(Mr -J^Tk^ 
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Figure 1. The pion scattering rate as a function of pion momentum with respect to 
medium calculated at T = 100 MeV and the highest estimates of chemical potentials 
allowed by data from SPS (left column) and RHIC (right column). Contributions to the 
total scattering rate from scattering on nucleons, anti-nucleons and pions are indicated. 
The lower row shows the baryonic and mesonic relative contributions. 

where the sum runs over all resonances in the ni scattering. Particle densities Pi(k) 
are assumed to be thermal. For realistic temperatures, only low-lying resonances are 
relevant and we include i = ir, N, N, K, p, A, A when summing in (|G|) over scattering 
partners. Since the freeze-out temperature is not known a priori, we scan three values: 
T = 90, 100, 120 MeV. Chemical potentials for pions are estimated from data on phase- 
space density, those for other species from measured ratios of dN/dy at mid-rapidity (see 
|| for more details on the estimates). 

Figure [I] shows the scattering rates as a function of the pion momentum relative to the 
heat bath calculated for SPS and RHIC, as well as the most important contributions to TZ. 
In spite of the strong increase of the pion phase-space density at RHIC |L|, the relative 
meson contribution to the total scattering rate does not grow proportionally since the 
7T7T cross-section is much smaller than the one for ttN. The total baryonic contribution 
changes very little since the smaller amount of baryons at RHIC is roughly compensated 
by anti-baryons. In summary, there is about 10% of the relative contribution shifted from 
baryons to mesons when going from SPS to RHIC. 

As seen in Figure 0, the scattering rate depends strongly on the assumed freeze-out 
temperature since the pion-nucleon scattering is typically dominated by total CMS energy 
below the A resonance peak. At higher temperature, average CMS energies move closer 
to the peak value thus leading to increased scattering cross-sections. Estimating the pion 
escape probability ([[]) from eq. (H) and £ = 0.08 fin -1 , one finds that a pion only has a 
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Figure 2. The pion scattering 
rate as a function of pion mo- 
mentum with respect to the 
medium calculated for differ- 
ent temperatures and a range 
of chemical potentials allowed 
by data (see New results 
of STAR collaboration cor- 
respond rather to the "high 
/x's" estimates. 



reasonable chance to escape (of order 50%) when the temperature is 100 MeV at most. 

We finally point out that the scattering rates in Figure |2] generically decreases with 
increasing pion momentum. This leads us to conjecture that high p±_ pions may decouple 
earlier, from hotter and smaller system. This may provide a novel as yet unexplored 
contribution to the observed M± dependence of HBT radii. 

The freeze-out criterion studied here differs from the Cooper-Frye formalism typically 
employed in hydrodynamical model studies. It will be interesting to explore to what extent 
this affects the predictions of hydrodynamical simulations for the collision evolution. 
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